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Abstract 

A preliminary measurement of the q 2 dependence of the D° — > K~e + v e decay rate is presented. 
This rate is proportional to the hadronic form factor squared, specified by a single parameter. This 
is either the mass in the simple pole ansatz m po i e = (1.854 ± 0.016 ± 0.020) GeV/c 2 or the scale 
in the modified pole ansatz a po ie = 0.43 ± 0.03 ± 0.04. The first error refers to the statistical, the 
second to the systematic uncertainty. 
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1 INTRODUCTION 



In exclusive semileptonic B decays, the accuracy of the determination of the contributing Cabibbo, 
Kobayashi and Maskawa (CKM) matrix elements |V^,| and \V u b\ is limited by the precision of the 
corresponding hadronic form factors participating in these decays. In c-hadron semileptonic decays, 
it can be assumed that \V CS \ and \ V c d\ are known (considering for instance that the Cabibbo matrix 
is unitary). Measurements of the hadronic form factors in such decays can thus be used to test 
predictions on their absolute value and q 2 (= (p e + p u ) 2 ) variation. In this expression, p e and p u 
are respectively the 4-momenta of the positron and of the neutrino. 

With the development of lattice QCD algorithms and the installation of large computing facil- 
ities, accurate values for hadronic form factors are expected from theory in the coming years. A 
comparison between these values and corresponding measurements of a similar and possibly better 
accuracy will validate these complex techniques and the uncertainties attached to the remaining 
approximations . 

In the present analysis the q 2 variation of the hadronic form factor, in the decay D° — > K~e + v e 5 , 
has been measured. Neglecting the electron mass, there is a contribution from a single form factor 
and the differential decay rate is a product of q 2 - and 8 e - dependent expressions where 9 e is the angle 
between the electron and the kaon in the ev e rest frame. Integrating over the angular distribution, 
the q 2 differential decay width reads: 



dr 

dq 2 



G 2 p 



^ = ^\Vcs\ 2 p 3 k (q z ) U(q z ) 



(i) 



As this decay is induced by a vector current generated by the c and s quarks, the q 2 variation 
of the form factor f+(q 2 ) is expected to be governed by the D* + pole. The following expressions 
have been proposed pQ: 



Mi 2 



and 12 



Uii 2 ) 



/+(Q) 

m' 2 , 

pole 



/+(0) 



(2) 



1 



7^ 



1 



(3) 



Equation (j2J is the "pole mass" and Equation © is the "modified pole mass" . Each distribution 
depends on a single parameter: m po i e and a po i e , respectively. In Equation @ rri£,* = (2.1121 ± 
0.0007) GeV/c 2 ) is the D* physical mass. In lattice QCD computations, usually a "lattice mass" 
value is used for mo* but the computed value for a po i e is expected to be directly comparable to 
the value extracted from the fit to data using expression ©. A recent result from lattice QCD 
computations 3^ gives dgjfc* = 0.50 ± 0.04. 

Radiative effects in D° — > K~e + v e decays have been simulated using the PHOTOS generator 
[1] and present measurements assume its validity. 



^Charge conjugate states are implied throughout this analysis. 
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2 THE BABAR DETECTOR AND DATASET 



Results given in this document have been obtained using BABAR data taken between February 2000 
and June 2002, corresponding to an integrated luminosity of 75 fb _1 and which comprises the Runl 
(18.6 fb _1 ) and Run2 (56.4 fb _1 ) data taking periods. In addition, Monte Carlo (MC) simulation 
samples of charm, 6-hadrons and light quark pairs, equivalent to about 1.3 times the data statistics 
have been used. 

The BABAR detector is described elsewhere [3] . Its most important capabilities for this analysis 
are charged-particle tracking and momentum measurement, charged tt/K separation, electron iden- 
tification, and photon reconstruction. Charged particle tracking is provided by a five-layer silicon 
vertex tracker and a 40-layer drift chamber. The Detector of Internally Reflected Cherenkov light 
(DIRC), a Cherenkov ring-imaging particle-identification system, is used to separate charged kaons 
and pions. Electrons are identified using the electromagnetic calorimeter, which comprises 6580 
thallium-doped Csl crystals. These systems are mounted inside a 1.5 T solenoidal superconducting 
magnet. 

3 ANALYSIS METHOD 

This analysis is based on the reconstruction of D* + mesons produced in continuum cc events and 
in which D* + — * D®ir + and the D° decays semileptonically. 

3.1 Candidate selection and background rejection 

Electron candidates are selected with a momentum larger than 0.5 GeV/c in the laboratory and, 
also, in the center of mass frame (cm). Muons are not used in this analysis as they have less purity 
than electrons and as we are not limited by statistics. The direction of the event thrust axis is 
taken in the interval | cos(# t h rU st)| < 0.6 to minimize the loss of particles in regions close to the 
beam axis and to ensure a good total energy reconstruction. 

Two variables, R2 (the ratio between the second and zeroth order Fox- Wolfram moments [H]) 
and the total charged and neutral multiplicity, have been used to reduce the contribution from B 
events based on the fact that the latter are more spherical. These variables have been combined 
linearly in a Fisher discriminant [Jj on which a selection requirement retains 71% of signal events 
and 10% of the B background. 

Charged and neutral particles are boosted to the center of mass system and the event thrust 
axis is determined. A plane perpendicular to the thrust axis, and containing the beam interaction 
point, is used to define two event's hemispheres. Each hemisphere is considered in turn to search 
for a candidate having an electron(±), a kaon(=p) and a pion(±) reconstructed in that hemisphere 
and with the relative charges as given within parentheses. A candidate triplet will be called also a 
right sign combination (RS), other different relative charge assignments will be called wrong sign 
(WS). 

To evaluate the neutrino (u e ) momentum, two constrained fits are used. In the first fit the 
(e + K~v e ) invariant mass is constrained to be equal to the D° mass whereas, in the second fit, the 
(e + K~v e TT + ) invariant mass is constrained to be equal to the D* + mass. In these fits, estimates 
of the D° direction and of the neutrino energy are included from measurements obtained from all 
particles registered in the event. The D° direction estimate is taken as the direction of the vector 
opposite to the momentum sum of all reconstructed particles but the kaon and the electron. The 
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neutrino energy is evaluated by subtracting from the hemisphere energy, the energy of reconstructed 
particles contained in that hemisphere. The energy of each hemisphere has been evaluated by 
considering that the total center of mass energy is distributed in two objects of mass corresponding 
to the measured hemisphere masses. The hemisphere mass is the mass of the system corresponding 
to the sum of the 4- vectors for particles contained in that hemisphere. Detector performance in the 
reconstruction of the D direction and missing energy have been measured using events in which 
the D° — > K~tt + and used in the mass constrained fits. 

After performing the first fit, the events with a x 2 probability larger than 10 -3 are kept. The 
D° 4-momentum is then obtained and the mass difference 5(m) = m(D°-K + ) — m(D°) is evaluated 
and is shown in Figure ^ 

Other particles, present in the hemisphere, which are not decay products of the D* + candidate 
are named "spectator" tracks. They originate from the beam interaction point and are emitted 
during hadronization of the created c and c quarks. The "leading" track is the spectator track 
having the largest momentum. Information from the spectator system is used, in addition to the 
one provided by variables related to the D° production and decay, to reduce the contribution from 
the combinatorial background. As charm hadrons take a large fraction of the charm quark energy, 
charm decay products have higher average energies than spectator particles. 

The following variables have been used: 

• the D° momentum after the first mass constrained fit; 

• the spectator system mass, which has lower values for signal events; 

• the direction of the spectator system momentum relative to the thrust axis; 

• the momentum of the leading spectator track; 

• the direction of the leading spectator track relative to the D° direction; 

• the direction of the leading spectator track relative to the thrust axis; 

• the direction of the lepton relative to the kaon direction, in the virtual W rest frame; 

• the charged lepton momentum in the cm frame. 

The first six variables depend on the properties of charm quark hadronization whereas the last two 
are related to decay characteristics of the signal. These variables have been linearly combined into 
a Fisher discriminant variable and events have been kept for values above 0, retaining 82% of signal 
events and rejecting 52% of background candidates. 

The remaining background from cc events can be divided into peaking and non-peaking cate- 
gories. Peaking events in 8(m) have a real D* + in which the slow tt + is included in the candidate 
track combination. The non-peaking background corresponds to candidates without a charged D* + 
slow pion. Using simulated events, these components are displayed in Figure ^ The MC values 
have been rescaled to the data luminosity, using the cross sections of the different components (1.3 
nb for cc, 0.525 nb for B + B~~ and B°B°, 2.09 nb for light uds quark events). 

To study the q 2 distribution, events with 5(m) below 0.16 GeV/c 2 and with a x 2 probability, 
of the second mass constrained fit greater than 1% have been selected. 
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Figure 1: 5{m) distributions from data and simulated events. The signal and the different back- 
ground components are indicated. MC events have been normalised to the data luminosity accord- 
ing to the different cross sections. The light quark component is very small and hardly visible. 



3.2 q 2 measurement 

The D° — » K~e + v e is a 3-body decay and its dynamics depends on two variables. In semileptonic 
decays these variables are usually taken as q 2 and cos(# e ). Throughout the present analysis, q 2 has 
been evaluated using q 2 = — pk) 2 {= (Pe + Pu e ) 2 ^j where pu and px are the four-momenta of 
the D and K meson, respectively. The differential decay rate versus these two variables factorises 
in two parts depending, respectively, on q 2 and 9 e . The variation of dT/dcos(6 e ) is fixed by angular 
momentum and helicity conservation and has a sin 2 9 e dependence. The aim of the present analysis 
is to measure the q 2 dependence of the form factor whose shape depends on the decay dynamics. 

After having applied the D° and D* + mass constrained fits, the values of the variables q 2 and 
cos(# e ) are obtained and will be noted as q 2 and cos(0 e ) r , respectively. Background contributions 
in these distributions have been evaluated from corresponding distributions obtained in simulated 
event samples, normalized to the same integrated luminosity analyzed in data. The measured q 2 
distribution is shown in Figure El The expected background level has been indicated and the fitted 
signal component, as obtained in the following, has been added. The q 2 distribution for signal 
events is obtained by subtracting the estimated background distribution from the total measured 
distribution. Using simulated events the reconstruction accuracy on q 2 has been obtained by 
comparing reconstructed (q 2 ) and true simulated (q 2 ) values. Figure El shows a double Gaussian fit 
to the q 2 resolution distribution. We obtain widths a± = 0.077 GeV 2 and oi = 0.276 GeV 2 for the 
two Gaussians, which are of similar area. 

The variation of the signal efficiency as a function of q 2 is shown in Figure 0J 
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Figure 2: Distributions of q 2 . Points with error bars are data. The distribution from background 
events is given by the blue histogram. The fitted signal component is overlayed. 



3.3 Unfolding procedure 

To obtain the unfolded q 2 distribution for signal events, corrected for resolution and acceptance 
effects, the Singular Value Decomposition (SVD) [Hj of the resolution matrix has been used in 
conjunction with a regularization which minimizes the curvature of the correction distribution. A 
SVD of a real m x n matrix A is its factorization of the form: 

A = U SV T , 

where U is an m x m orthogonal matrix, V is an n X n orthogonal matrix, while S is an m x n 
diagonal matrix with non- negative elements: 

Sii = Si ^ 0. 

The quantities Si are called Singular Values (SV) of the matrix A. If the matrix of a linear system 
is known with some level of uncertainty, and some SV of the matrix are significantly smaller than 
others, the system may be difficult to solve even if the matrix has full rank, and SVD suggests a 
method of treating such problems. We will assume that the SV values Sj form a non-increasing 
sequence. 

The method uses binned distributions. It needs, as input, a 2-d array which indicates how 
events generated in a bin in q 2 are distributed over several bins in q 2 . To be able to correct for 
acceptance effects one needs also the initial distribution in q 2 , as given by the generator. The 
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Figure 3: Distribution of the difference between the real and the reconstructed value of the q 2 
variable (Aq 2 = q 2 — q 2 ) as obtained from simulated events. The fitted curve is the sum of two 
Gaussian distributions. 



number of events estimated in each bin of the q 2 variable ( n\ 



estimated\ i j 

2 can be expressed as: 



Ei j i,estimated 

7=1 



where w 3 is an unknown deviation of the number of simulated events in bin j of the q 2 variable, 
from the initial Monte Carlo value. Aij is the actual number of events which were generated in 
bin j and ended up in bin i. The elements of the vector are determined by minimizing a x 2 
expression obtained by comparing estimated and measured number of events in each q 2 bin: 



E 



n" 



mm, 



where a n i is the corresponding uncertainty on the measured number of events in bin i. The 



minimization leads to the system: 



Aw 



n q 2 



(4) 



where A 



Kl 



Aij/a n i and h q 2 = n l 2 ja nl 

Qr It 

The exact solution of Equation @ will again most certainly lead to rapidly oscillating distribu- 
tion. This spurious oscillatory component should be suppressed and this can be achieved by adding 
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Figure 4: q 2 dependent efficiency after all selection requirements, from (Runl) simulated signal 
events. 



a regularization term to the expression to be minimized: 

(^Aw — nj (^Aw — nj + r (Cw) T Cw = min. 

Here C is a matrix which suppress solutions w having large curvatures and r determines the relative 
weight of this condition. The curvature of the discrete distribution Wj is defined as the sum of the 
squares of its second derivatives: 

i( w i+i ~ w i) ~ ( w i ~ w i-i)? ■ 

i 

The unfolded q 2 distribution is given by the bin-by-bin product of w and of the initial Monte 
Carlo distribution. This approach provides the full covariance matrix for the bin contents of the 
unfolded distribution. Its inverse should be used for any x 2 calculation involving the unfolded 
distribution. 

For a square response matrix, Singular Values (SV) correspond essentially to eigenvalues of 
this matrix. The number of significant SV has been obtained by looking at their distribution in 
experiments generated using a toy simulation of the experimental conditions, and keeping those 
which are above a plateau |Hj. According to this procedure, four values can be retained in the 
present analysis. 

Unfolded q 2 distributions, obtained by dividing the total sample into four similar data sets 
(each corresponding to about 21k events) have been analyzed. These distributions, and the total 
resulting sample, have been fitted using the pole and the modified pole models. Corresponding 
values for the parameters are given in Tabled 
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It has to be noted that obtained values for m po i e and a po i e are independent of a choice for 
a number of SV. In addition selecting a fixed number of SV results in introducing biases on the 
values evaluated in each bin for the unfolded distribution. Thus, at present, when providing the q 2 
dependence of the hadronic form factor and the corresponding error matrix we have kept all SV 
(see also Section I4.5|l . 

Table 1: Fitted values for the parameters corresponding respectively to a pole mass and a modified 
pole mass model for the form factor, obtained in different samples of similar statistics and for the 
total. The last column indicates the fraction of background events. 



Data sample 




(GeV/c 2 ) 


O^pole 


bckg. fraction 


Runl 


1.80< 


3± 0.027 


0.54 ±0.06 


16.3% 


Run2-1 


1.88« 


3 ± 0.033 


0.37 ±0.07 


16.6% 


Run2-2 


1.847 ±0.032 


0.43 ± 0.07 


16.4% 


Run2-3 


1.864 ±0.033 


0.41 ±0.07 


16.7% 


All 


1.854 ±0.016 


0.43 ± 0.03 


16.5% 



4 SYSTEMATIC STUDIES 

Systematic uncertainties originate from non-perfect simulation of the charm fragmentation process 
and of the detector response, uncertainties in the control of the background level and composition, 
and the unfolding procedure. Differences between data and Monte Carlo in quantities used in the 
analysis may result in biases that need to be corrected. 

These effects have been evaluated by redoing the unfolding procedure after having modified the 
conditions that correspond to a given parameter entering in the systematics and evaluating the 
variation on the value of the fitted m po i e and a po i e parameters. 

For some of these studies dedicated event samples have been used in which a D° is reconstructed 
in the K~ir + or K~-k + it® decay channel. 

4.1 Systematics related to c-quark hadronization 

The signal selection is based on variables related to c-quark fragmentation and decay properties of 
signal events. As we measured differences between the hadronization properties of c-quarks events 
in real and simulated events, a g 2 -dependent difference in efficiency between these two samples is 
expected. To have agreement between the distributions measured in real and simulated events, for 
the different variables entering into the Fisher analysis a weighting procedure has been used. These 
weights have been obtained using events with a reconstructed D° decaying into K~ir + or K~it + it® . 
The data-MC agreement of the measured distributions indicates that uncertainties related to the 
tuning are below 5 MeV/c 2 on m po i . These corrections correspond to a 16 MeV/c 2 displacement 
on m po i e . 
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4.2 Systematics related to event selection and analysis algorithm 

Another study has been done by analyzing D° — > K^ir + 7r° events as if they were K~e + v events. 
The two photons from the 7r° are removed from the particle lists and events are reconstructed using 
the algorithm applied to the semileptonic D° decay channel. The "missing" tt° and the charged 
pion play, respectively, the role of the neutrino and of the electron. In the two mass constrained fits, 
to preserve the correct kinematic limits, it is necessary to consider that the "fake" neutrino has the 
7T° mass and that the "fake" electron has the 7r + mass. With these events one can measure possible 
g 2 -dependent effects related to some of the selection criteria and the accuracy of the measurement 
of q 2 . 

In these studies the "exact" value of q 2 is defined as q 2 = (p n + +p n o) 2 . Results obtained 
using real and simulated events are compared to identify possible differences. Fractions of real and 
simulated events are selected using the Fisher discriminant variables defined to reject background 
events, the % 2 probability from the two fits, the 5(m) selection, and the requirement for having at 
least one spectator particle. In this selection, the Fisher variable used against the cc background 
does not include cos(# e ) and p e . This is because the angular distribution in K~ir + ir° events is close 

to COS 2 (#„-). 

There is no evidence for a statistically significant different behavior in data as compared with 
the simulation. The statistical accuracy of this comparison, based at present on the Runl data 
sample, corresponds to an uncertainty on the pole mass of ±15 MeV/c 2 . 

4.3 Detector related systematics 

To search for differences between real and simulated events, distributions of Aq 2 = q 2 — q 2 , obtained 
by selecting D° — ► K~ir + ir° events in a given bin of q 2 have been compared. One notes that these 
distributions are systematically slightly narrower for simulated events. We have evaluated this effect 
to be lower than 5% and have redone the exercise after adding a smearing on q 2 , in simulated events. 
We obtain a +4 MeV/c 2 increase on the fitted pole mass and assign this value as a systematic error. 

Effects from a momentum dependent difference between real and simulated events on the 
charged lepton efficiency reconstruction and on the kaon identification have been evaluated. Such 
differences have been measured for electrons and kaons using dedicated data samples. When ap- 
plied, the corresponding corrections induce an increase of the fitted pole mass by +2 MeV/c 2 and 
+4 MeV/c 2 , for electrons and kaons respectively. A relative uncertainty of 30% has been assumed 
on this correction. 

4.4 Background related systematics 

The background under the D* + signal has two components which have respectively a peaking and 
a non-peaking behavior. 

The peaking background arises from events with a real D* + in which the slow ir + is included in 
the candidate track combination. Its main components, as expected from the simulation, comprise 
events with real or fake K~ or positron. They mainly correspond to D° — > K~iv e + iy e and D° — > 
K~it + it events where, for the latter, the 7r° has a Dalitz decay or a converted photon. 

The non-peaking background originates from non-cc events and from continuum charm events 
in which the tt candidate does not come from a decaying D* + . RS combinations, for 5(m) > 
0.18 GeV/c 2 and WS combinations can be used to compare data and simulated event rates. It has 
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been verified that 5{m) distributions were in agreement in the tail of RS events and reproduce well 
the level of WS combinations (see Figures ^ and 0) . 
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Figure 5: <5(m) distributions from data (points with error bars) and simulated events (shaded 
histogram) in WS K~e + ir~ events. 



Table 2: Systematics corresponding to uncertainties on the level of the different background sources. 



Source 


variation 


5 (m P oi e ) 


<5 ("pole) 




(%) 


(MeV/c 2 ) 






±9 


±4.4 


=F0.0087 


other peaking 


±10 


±1.2 


T0.0023 


cc non-peaking 


±10 


±3.0 


=F0.0062 


B+B- 


±10 


±2.6 


^0.0050 


B°B° 


±10 


±4.3 


^0.0084 


uds 


±10 


±0.5 


=F0.0010 


total 




±7.4 


±0.0147 



As outlined in Tabled each background component has been varied by ±10% ( a value which 
is larger than observed differences in the analyzed distributions), apart for the D° — > K~i:^e r v e 
component for which we use ±9%, corresponding to the combination of the PDG uncertainty with 
a recent measurement of this decay channel by the CLEO-c experiment j^j. 
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4.5 Control of the statistical accuracy and of systematics in the SVD approach 

For a fixed number of singular values, we verify that the statistical precision obtained for each 
binned unfolded value is reasonable and that biases generated by having removed information are 
under control. 

This has been achieved by fitting the same model for the form factor, directly on the measured 
distribution, after having included resolution effects, and then by comparing the uncertainty ob- 
tained on the extracted parameter with the one determined by fitting the unfolded distribution 
using that model. These studies have been complemented by toy simulations. 

One observes that the error obtained from a fit of the unfolded distribution is underestimated 
by a factor of ~ 1.3, which depends on the statistics of simulated events. Pull distributions indicate 
also that unfolded values in each bin are biased. The importance of the biases decreases as the 
number of kept singular values increases. When the number of SV is equal to the number of bins no 
biases are present. When the number of SV is smaller, these effects are corrected using pull results. 
It has been verified that these corrections, obtained for a given value of the pole mass, do not bias 
results expected in experiments generated with a different value for m po i e . Corrections defined with 
"T-poie = 1-8 GeV/c 2 induce a shift below 4 MeV/c 2 when applied on experiments generated with 
m po ie = 1.7 or 1.9 GeV/c 2 . 

4.6 Systematics summary 

The systematic uncertainties evaluated on the determination of m po i e and a po i e are summarized in 
Table H 

Table 3: Summary of corrections and systematics on the fitted parameters. The sign of the correc- 
tion is given by the sign of the difference between the fitted values after and before applying the 
correction. 



Source 


8 ("Vie) 






(MeV/c 2 ) 




c-hadronization tuning 


-16 ±5 


±0.03 ±0.01 


reconstruction algorithm 


±15 


±0.03 


resolution on q 2 


+4 ±4 


-0.01 ±0.01 


particle ID 


+6 ±2 


-0.02 ±0.01 


background control 


±7 


±0.02 


unfolding method 


±5 


±0.01 


total 


-6 ±20 


±0.01 ±0.04 



The systematic error matrix for the ten unfolded values has been also computed by considering, 
in turn, each source of uncertainty and by measuring the variation, Si, of the corresponding unfolded 
value in each bin (i). The elements of the error matrix are the sum, over all sources of systematics, 
for the quantities 5{ ■ 5j. The total error matrix has been evaluated as the sum of the matrices 
corresponding respectively to statistical and systematic uncertainties. 
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5 RESULTS 



In the present measurement the fit to a model has been done by comparing the number of events 
measured in a given bin of q 2 with the expectation from the exact analytic integration of the 
expression p\{q 2 ) |/+(g 2 )| over the bin range. The normalisation is left free to float and the 
weight matrix evaluated from SVD is used. From the measured integrated decay spectrum in each 
bin, the value of the form factor |/ + ((/ 2 )| has been evaluated at the bin center. By convention it 
has been assumed that |/+(0)| = 1. 

The unfolded q 2 distribution, for signal events and keeping four SV, is given in Figure El where 
it has been compared with the fitted distributions obtained using the pole and the modified pole 
ansatz. 
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Figure 6: Comparison between the normalised unfolded q 2 distribution of the form factor obtained 
in the present analysis, keeping four SV, and those corresponding to the two fitted models. Lower 
plots give the difference between measured and fitted distributions. Vertical lines correspond to 
statistical uncertainties. 



Results for the values of m po i e and a po ie are independent of the choice of a given number of 
kept SV. The number of SV corresponds to the number of independent combinations of the ten 
bins used in the analysis. To minimize correlations between bin's content and to compare with 
the q 2 distributions obtained in other experiments the distribution obtained using ten SV has been 
provided. The content of each bin is given in Table 0] The corresponding statistical and total error 
matrices are given also in that table and the q 2 variation of the hadronic form factor is displayed 
in Figure 

Similar measurements of the q 2 dependence of |/+(g 2 )| have been obtained by several experi- 
ments; recent published results are summarized in the following. 
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Table 4: Unfolded distribution using ten SV. The second line of this table gives the integrated 
values of the differential decay branching fraction over 0.2(GeV 2 ) intervals (quoted on the first 
line). The total distribution has been normalised to unity for q varying from to 2 (GeV 2 ). The 
statistical error matrix corresponding to the above measurements is then provided. On the diagonal 
is given the uncertainty on each measured value. Off-diagonal terms correspond to the correlation 
coefficients. Including systematics, the total error matrix is then provided. 



2 u- 

q bin 


[0, 0.2] 


[0.2, 0.4] 


[0.4, 0.6] 


[0.6, 0.8] 


[0.8, 1.0] 


[1.0, 1.2] 


[1.2, 1.4] 


[1.4, 1.6] 


[1.6, 1.8] 


[1.8, 2.0] 


fraction 


0.1999 


0.1791 


0.1606 


0.1489 


0.1007 


0.0962 


0.0570 


0.0363 


0.0201 


0.0012 


stat. 


0.0040 


-.63 


.24 


-.094 


.034 


-.0065 


-.0022 


0.0006 


-.0025 


.0002 


error 




0.0072 


-.66 


.27 


-.10 


.025 


.0015 


-.0003 


.0008 


-.0001 


and 






0.0090 


-.69 


.29 


-.081 


.011 


-.0074 


.0022 


.0004 


correl. 








0.0081 


-.67 


.24 


-.060 


.025 


-.0063 


.0019 












0.0071 


-.64 


.23 


-.083 


.021 


-.0066 














0.0056 


-.64 


.25 


-.064 


-.0035 
















0.0044 


-.60 


.15 


.022 


















0.0032 


-.44 


-.050 




















0.0018 


-.059 






















0.00054 


total 


0.0045 


-.48 


.23 


-.097 


-.0084 


-.073 


-.079 


-.082 


-.12 


-.044 


error 




0.0073 


-.65 


.27 


-.11 


.0038 


-.023 


-.027 


-.039 


-.016 


and 






0.0090 


-.69 


.28 


-.086 


.0027 


-.016 


-.013 


-.0071 


correl. 








0.0081 


-.67 


.24 


-.055 


.029 


-.0048 


.0006 












0.0071 


-.62 


.24 


-.067 


.039 


-.0003 














0.0057 


-.60 


.27 


-.028 


.0075 
















0.0045 


-.55 


.19 


.041 


















0.0032 


-.36 


-.025 




















0.0019 


-.0006 






















0.00055 
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Data taken at the T(4S) energy and corresponding to an integrated luminosity of 7 fb _1 have 
been analyzed by CLEO III [111 . Only three q 2 bins have been used because of the rather poor 
resolution obtained on this variable. 

The FOCUS fixed target photo-production experiment has collected 12840 signal events for the 
decay chain: D*+ -> D°tt + ; D° -> R-fi+u^ 12 . The q 2 resolution has a r.m.s. of 0.22 (GeV/c) 2 . 

The BELLE i?-factory experiment has done a measurement by reconstructing all particles in 
the event, but the neutrino, which parameters can be obtained using kinematic constraints. As a 
result, they achieve a very good reconstruction accuracy of about 0.010 (GeV 2 ) on q 2 but at the 
price of a low overall efficiency. Analyzing 282 ft -1 integrated luminosity, they select about 2500 
> + f M and D° — ► K e + v e events with low background levels. 

In Table |S] results obtained by the different collaborations have been summarized when fitting 
a pole mass and a modified pole mass q 2 distributions for the form factor. 



Table 5: Fitted values for m po i e and a po i e models for the form factor. 



Experiment 


m poie (GeV/c 2 ) a po ie Statistics 


CLEO III UU 


1.89 ±0.05t£^ 0.36±0.10t&$ 7 ft.- 1 


FOCUS | 


12 


1.93 ± 0.05 ± 0.03 0.28 ± 0.08 ± 0.07 13k events 


BELLE 




1.82 ± 0.04 ± 0.03 0.52 ± 0.08 ± 0.06 282 ft" 1 


BABAR 


1.854 ±0.016 ±0.020 0.43 ± 0.03 ± 0.04 75 ft" 1 



Results obtained in this analysis on the q 2 variation of the hadronic form factor have been 
compared in Figure with FOCUS [12] measurements and with the lattice QCD [2] evaluation. 

The present result is valid in the limit that in real events the radiative component has the same 
characteristics as in the simulation which uses the PHOTOS generator program. With the level of 
accuracy of present measurements better studies of radiative effects are needed. 

6 SUMMARY 

Fitting the pole mass and the modified pole mass ansatz to the measurements, we obtain preliminary 
values for the single parameter that governs their q 2 dependence: 

mpol0 = (1.854 ± 0.016 ± 0.020) GeV/c 2 (5) 
a pole = 0.43 ± 0.03 ± 0.04 

where the first error is statistical and the second systematic. The effective m po i e value is rather 
different from m n *+ indicating large corrections to naive expectations. In the modified pole model 
this can be interpreted as evidence for the contribution from other vector states of invariant mass 
higher than the mass. 

The value measured for a po i e agrees, within errors, with the one obtained from lattice QCD pJJ: 
a pole Ce = O-^O i 0.04. We provide also a preliminary q 2 distribution of the form factor, corrected 
for effects from reconstruction efficiency and finite resolution measurements. 
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Figure 7: Comparison between the measured variation of < \f+(q 2 )\ > obtained in the present 
analysis and in the FOCUS experiment. The band corresponds to lattice QCD (3] expectations. 
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